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Upper Limits on the Pulsed VHE 7-ray Emission from 
Two Young Pulsars Investigated with the High Energy 
Stereoscopic System 
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Abstract. The High Energy Stereoscopic System 
(H.E.S.S.) is a system of four, imaging, atmospheric 
Cherenkov telescopes in Namibia, designed to detect 
very-high-cnergy 7 rays above ~ 100 GcV. During 
2002-2003, H.E.S.S. collected data from two, young and 
energetic radio pulsars: the Crab and PSR B1706— 44. 
We searched for pulsations at the lowest energies that 
H.E.S.S. is capable of detecting, aiming at a detection 
that would potentially differentiate between the two 
popular models of pulsar high-energy emission: the Polar 
Cap and the Outer Gap. No evidence for pulsed emission 
was found in the data, and upper limits were derived to 
a 99.95% confidence level. Our assumptions and upper 
limit values for the two pulsars are reported. 



1. Introduction 

The Crab pulsar (PSR B0531+21) is a well-known source 
of electromagnetic pulsed emission in a broad range of fre- 
quencies, from radio to high-energy (HE > 100 MeV) 7 
rays (Thompson 1999). Virtually all efforts to detect this 
emission has led to success, and this pulsar has been justi- 
fiably characterized as the "standard candle" . EGRET has 
detected its pulsed signature up to ~ 20 GeV (Nolan et al. 
1993). Yet, Crab's pulses become undetectable above this 
range, and so far there has not been a verifiable detection 
in very high energies (VHE > 100 GeV). 

PSR B1706— 44 is a young, southern-hemisphere pul- 
sar that lies relatively close to the Galactic centre. It was 
first discovered in radio, but it also appears as a strong 
emitter in HE 7 rays (Johnston et al. 1992; Thompson et 
al. 1996). The detection with EGRET in the latter energy 
range has provided the only confident high-energy profile 
so far, although a strong indication of pulsed emission, to 
a 4(j level, was also found in X-ray data from Chandra 
(Gotthelf et al. 2002). Above ~ 20 GcV, only upper limits 
exist for this pulsar. 

A few of the above pulsars' properties are summarised 
in TableHJ 
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Fig. 1. Geometrical model of a pulsar's magnctosphcrc 
with the locations of emission according to the Polar Cap 
and Outer Gap models. 



2. Theory 

The absence of pulsed emission above EGRET 's range is 
consistent with two popular classes of theoretical mod- 
els, the Polar Cap (PC) and the Outer Gap (OG). They 
predict that the high-energy spectra of pulsars should di- 
minish rapidly beyond a certain cut-off energy. 



2.1. Polar Cap models 

According to PC models, the observed pulsed emission is 
generated close to the pulsar's surface, above the polar 
caps (Fig. [I}. High-energy 7-ray production is triggered 
mainly via synchro-curvature radiation from relativistic 
e~ and e + . Such models predict HE spectra with steep cut- 
offs at a few GeV (Sturrock 1971; Ruderman & Sutherland 
1975). 
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Table 1. The properties of the Crab and PSR B1706-44 
pulsars 





Crab 


PSR B 1706 44 


P(s) 


0.033 


0.102 


P (x 1(T 15 s s" 1 ) 


421 


93 


log T (y) 


3.1 a 


4.2 


logBs (G) 


12.6 


12.5 


E (x 10 36 erg s" 1 ) 


450 


3.4 


e 7 % (> 100 MeV) b 


~ 0.01 


~ 0.1 



a This is the dynamic age, |(P/P); Crab's true age is 952 y. 
b HE 7-ray conversion efficiency: L 7 (> 100 MeV)/P 



2.2. Outer Gap models 

The OG models, on the other hand, place the site of emis- 
sion in the outer magnctosphcrc, inside charge-depleted 
regions (gaps) near the null-charge surface (Fig. |T|). The 
main radiation mechanism is synchro-curvature radiation 
from e~-e + pairs accelerated along the magnetic field lines 
inside the gaps. Furthermore, OG models predict an addi- 
tional inverse Compton component from the upscattering 
of low-energy photons (i.e. IR, optical, X-ray) by acceler- 
ated e~-e + pairs. In general, OG models produce gentler 
spectral cut-offs than those of PC models (Cheng, Ho & 
Ruderman 1986; Chiang & Romani 1992). 



3. Observations 

H.E.S.S. is a system of four imaging Cherenkov telescopes 
at the Khomas highland of Namibia (Aharonian et al. 
1997b). The array is capable of detecting, at 5er signifi- 
cance, TeV point sources that are 100-times weaker than 
the Crab nebula, in 50 h exposure time. At its threshold 
(~ 100 GeV), H.E.S.S. can still detect sources that are 
10-times weaker than the Crab nebula, given the same 
exposure time. 

In 2003, H.E.S.S. performed stereoscopic observations 
of the Crab, using three of the four telescopes of the array. 
During that time, the source was observable between 45° 
and 60° zenith angle (Z.A.). For our analysis, we selected 
only data fulfilling quality criteria, like good weather con- 
ditions and high trigger rate (R > 150 Hz). The total 
exposure time of the selected data set amounted to T ~ 
4.5 h. 

In 2002, H.E.S.S. performed single-telescope observa- 
tions of PSR B1706— 44, using the first telescope available 
then. The Z.A. of the observations was < 30°. Based on 
the same quality criteria, an extensive data set comprising 
28 h of data was selected for analysis. 



4. Data analysis 

4-.1. Low-energy cuts 

Our analysis targeted the lowest energies that H.E.S.S. 
could detect in the above observations, in an effort to in- 
tercept the tail-end of the pulsed emission observed with 
EGRET. For that purpose, we applied specially tailored 
event-selection cuts to the 7-ray images: we rejected all 
events whose total photo-electron content was above 100 
and accepted only events which laid closer than 18 mrad 
from the source position. After the application of cuts, the 
background level was ~ 10 6 events for both data sets. 

4-2. Temporal analysis 

Using standard epoch- folding techniques (EF), the times 
of arrival (TOAs) of the events were folded into phases. For 
that purpose, an accurate, contemporaneous ephemeris 
was needed: for the Crab, it was obtained from Jodrcll 
Bank, and for PSR B1706-44 we used the GRO ephemeris 
of the Australian Pulsar Timing Archive. The resulting 
phasograms after binning the event phases across two pul- 
sar periods are shown in Fig. [2] Both phasograms ap- 
pear 'flat' and consistent, within the errors shown, with 
Poissonian-background fluctuations. 

Despite the apparent uniformity, we decided to search 
for a potential signal, hidden in the large background, with 
more sensitive tests for uniformity: those included e.g. the 
-ff-test and the x 2 -test (de Jager et al. 1989). The chance 
probabilities from the H-test are shown in Table HJ On 
the whole, all chance probabilities from all our tests were 
consistent with the "null hypothesis" : i.e. the lack of pe- 
riodicity in the data. 

4-3. Upper limits on the pulsed flux 

The uniformity tests did not reveal significant deviations 
from the background. However, this could be due to the 
relatively weak signal compared to the background fluc- 
tuations. Hence, knowing the magnitude of the latter, one 
can set upper limits (hereafter ULs) to the value of the 
pulsed flux. 

We decided to calculate 3er (or 99.95%) ULs on the 
7-ray counts, using the method of O. Helene (1983). He- 
lene's method requires prior knowledge of the phase ranges 
where pulsed emission occurs ("peak areas"), across the 
pulsar's period. For the Crab and PSR B1706-44, the 
peak areas are known from HE observations with EGRET 
(see Fig. [2]); we assumed that those ranges remain the 
same in the VHE regime. 

In order to translate the above ULs to flux ULs, we 
utilised Monte Carlo simulations to estimate the tele- 
scope's effective area (A c s) for those events passing our 
cuts. The A c g(E) functions were derived for Z.A. corre- 
sponding to the observed data. 
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Fig. 2. Phasograms of the folded events from 4 h of Crab 
pulsar data and 28 h of PSR B1706-44 data. The shaded 
areas correspond to the phase ranges where significant 
pulsed emission was seen with EGRET. 



4- 4- Energy thresholds 

An important parameter that would help quantify how low 
in energy our cuts have allowed us to probe is the energy 
threshold of the data sets after cuts. Typically, the energy 
threshold for a particular sample of events — with a cer- 
tain range of energies, incident angles, distances from the 
telescope, etc. — represents the energy below which the 
detector's A e g for these events declines rapidly. More ac- 
curately, the energy threshold is defined as the maximum 
of the differential trigger rate: the latter being equal to the 
product between A e g and the differential flux spectrum, 
dN/dE. For the latter, we simply assumed an extrapola- 
tion of EGRET's HE power laws, dN/dE oc E~ v , beyond 
~ 20 GcV. The assumed spectral indices, v, are shown in 
Table El 

By convolving A e g with dN/dE, we plotted the dif- 
ferential trigger rates, dR/dE, for the two pulsars (see 
Fig.[3|). It can be seen that, although there is a clear max- 
imum at an energy E t h in both cases, the low-energy tail 
below -Eth accounts for a large fraction of the total trig- 
ger rate; and, clearly, this part must be considered in the 
calculation of our ULs. Hence, we decided to use a more 
representative definition of E t h in our analysis: we set E' th 
to be the energy at which dR/dE drops to 1/e of the max- 
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Fig. 3. Differential trigger rate functions (dR/dE) for the 
Crab and PSR B1706— 44, as a function of photon energy, 
for the assumed power-law spectra with v = 2.08 and v = 
2.1, respectively. Also shown are the standard definition 
of the energy threshold (-E'th), as the maximum of dR/dE, 
and our definition of E t ' h , at 1/e of the maximum dR/dE. 



Table 2. Analysis parameters and resulting fluxes for the 
Crab and PSR B1706-44 





Crab 


PSR B 1706 44 


T(h) 


4.5 


28 


Z.A. 


45-60° 


< 30° 


-^background (Hz) 


98 


45 




0.29 


0.008 


s (%py 


10 


30 




2.08 


2.1 


E' th (GeV) 


232 ± 51 


75 ±12 


A eS (K b ) (m 2 ) 


10 4 


400 


F u i(> E' th ) (x 10- 10 crn^s- 1 ) 


1.8 


25 



a Background trigger rate after cuts 
b Chance probability under the //-test 

c Duty cycle of the pulsed emission observed with EGRET 
above 100 MeV 



imum value, as one moves down in energy. This led to the 
thresholds shown in Table [2j 



5. Results and conclusions 

Our ULs for the Crab and PSR B 1706-44 are shown 
in Fig. |4] and Up to ~ 20 GeV, there exist confident 
detections with EGRET, whereas at higher energies the 
emission has only been constrained with ULs from various 
VHE ground-based experiments. It can be seen that our 
ULs are well-below the energies of those after the standard 
H.E.S.S. analysis and below many of other experiments. 
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Fig. 4. The derived upper limit (UL) on the differential 
pulsed 7-ray flux of the Crab pulsar, after the applica- 
tion of our low-energy cuts ("This work"). Flux ULs from 
various experiments are shown in light grey, with differ- 
ent symbols (Lessard et al. 2000; Dowthwaite et al. 1984; 
Tanimori et al. 1998; Aharonian et al. 2004; Lopez 2005). 
The UL at 60 GeV has been derived from observations 
with CELESTE (de Naurois et al. 2002). All ULs la- 
belled "H.E.S.S." have been derived from the standard 
H.E.S.S. analysis (solid black squares) (Schmidt 2005). 
Also, the Polar Cap (PC) and Outer Gap (OG) model 
spectra for this pulsar arc plotted with a smooth, light- 
grey and a stepped, grey line, respectively (Schmidt 2005; 
Hirotani & Shibata 2001). The squigglcd-dottcd line is a 
power-law fit to EGRET 's data points, extrapolated to 
very high energies (Nolan et al. 1993). 



Unfortunately, due to the large background inherent in 
our low-energy analysis, the corresponding ULs are 2-3 
orders of magnitude larger than those at higher energies. 

In the case of the Crab pulsar, we can conclude that 
the derived UL excludes confidently the possibility of a 
single power law from EGRET 's range to at least ~ 230 
GeV, thus verifying the indications for a cut-off already 
observed at the top energy bin of EGRET. However, for 
PSR B1706— 44, the large UL value prevents a similar con- 
clusion. 

Alongside the experimental results, we have also plot- 
ted the spectral predictions from a PC and OG scenario 
for these two pulsars. One of the main differences between 
PC and OG spectra is the much steeper cut-offs of the for- 
mer compared to the latter; this could provide the key to 
determining which model describes pulsar emission best. 
However, it is clear that our ULs, although nearer to the 
spectral cut-offs than previous H.E.S.S. results, still lack 
the sensitivity in terms of flux at the required low energies; 
hence the problem stands. 



> 
g 

UJ 

TJ 
Z 



10" 
10" 
10" 
10" 



10" 



10" 



10" 
10" 



10 



■10 



10" 



10 







EGRET 




This work Durham 






vT^g* , H.E.S.S. 


■ 


--Q.1 
* ***v*v * E 

PC 






OG 








PSR B 1706-44 








) 2 10 3 10 4 




10 6 10 6 


10 7 



E (Me VI 



Fig. 5. The derived upper limit (UL) on the differential 
pulsed 7-ray flux of PSR B1706— 44, after the application 
of our low-energy cuts ("This work"). The UL from ob- 
servations with Durham's Mark 6 telescope is shown with 
a grey, empty circle (Chadwick ct al. 1998). All ULs la- 
belled "H.E.S.S." have been derived from the standard 
H.E.S.S. analysis (solid black squares) (Schmidt 2005). 
Also, the Polar Cap (PC) and Outer Gap (OG) model 
spectra for this pulsar are plotted with a smooth, light- 
grey and a stepped, grey line, respectively (Schmidt 2005; 
Hirotani 2003). The squigglcd-dottcd line is a power-law 
fit to EGRET's data points, extrapolated to very high 
energies (Thompson et al. 1996). 



Future experiments, like H.E.S.S. Phase II, MAGIC 
and GLAST, will almost certainly fill the energy gap be- 
tween HE and VHE observations and, consequently, piece 
the puzzle of pulsar 7-ray emission. 
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